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A Mathematical Expression for the Growth 
of Trees in their Dependence on Time and 


Density of Stocking 


Afforestation in South Africa is mainly concerned with 
evenvaged coniferous plantations. In this paper an attempt is 
made to derive a mathematical expression for the basal area 
increment in such stands in relation to the density of stocking 
and the time. 


The growth of a tree is governed in the main by the 
following four factors — 


(a) The tree itself, ie., mainly the species. 
(b) The soil guality. 

(c) The water supply. 

(d) The climatic conditions. 


If we could ensure these factors to be the same the 
probability is that all trees would grow alike. 'The diferen- 
tiation in growth that obtains in practice may be ascribed to 
a change in some or all of these factors. A tree will have its 
optimum growth when these factors are operating in a way 
determined by the demands made by the nature of the treer 
species. When this balance is departed from the optimum 
rate of growth is affected and we may then speak of restricted 
growth. 'This balance may be disturbed by natural conditions, 
e.g. prolonged drought, or by mutual competition. 'This 
competition again is a function of the density of stocking. It 
becomes practically infinite and tends to stifle all growth when 
the seedlings are contiguous. (On the other hand when trees 
are so far apart that they do not influence one another's growth, 
they are termed free-growing. 


We may define the growing space A of a tree as the 
region from which it abstracts its food-supply. 'The growing 
space of the freesgrowing tree varies with the species, the soil, 
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the climate, the watersupply and its age. lt may sometimes 
cover a comparatively large fraction of an acre. 


For the purpose of this investigation we will disregard 
annual fluctuations in rainfall and climate and hence assume 
the water-supply for a tree on a growing space A as constant 
from year to year. ln so far as the watersupply is constant, 
the food-supply will also be so, since the tree's growth is 
arrested by a failing of the limiting factor, water. 'The water 
supply per unit of time (one year) is proportional to A, and 
the food-supply per annum available to the tree is thus b'A, 
where b! is a constant of proportionality. If b'A is less than 
the maximum redguired for the same tree when freegrowing 
in the same locality, all this food is abstracted in the course of 
one year and enters into the metabolism of the tree, whereby a 
certain guantity is built up into the body of the tree, roots, stem, 
bark and leaves. 'The rest of the food abstracted is discharged by 
the tree in the form of transpiration and the dropping of bark, 
leaves and twigs. 'The discharge we assume proportional to 
the crown of the tree, ie., the number of leaves of the tree. 
If we assume with Matthews (3) that by weight the leaves 
of a tree bear a definite linear relationship to the stem and 
branches of the tree, the food-supply discharged will be 
proportional to the total weight of the stem and branches, 
which again may be assumed to be proportional to its volume. 
The food discharged is thus fv, where v is the volume of the 
stem and branches, and f a factor of proportionality. 'The 
same relation will hold, with a slightly different f, if we neglect 
the weight of the branches, compared to the stem, which 
practically obtains in a forest. 'The food-supply used up and 
AE up into wood, bark, roots and leaves per unit of time 
1$ thus 


DA — fo 


If a volume dv of wood is thus put on in the time dt we must 
have— P 


dy co (HA — fo) di 
—k (BA — fo) di 


where k is a factor of proportionality. 
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We want to reduce this eguation to guantities which are 
more accessible to measurement than volumes. We are enabled 
to do so by assuming with Kopezky (2) that the volumes 
of trees are in a linear relation to their basal area— 


ie. v is proportional to a 


— ca Say 
dv — cda 


On substitution above we get— 
da —E (BA — fea) d 


AG AE HA —f) 


DAE AE Ad 


This is a differential eguation for the growth of trees in a 


- da 
forest, in which we may regard : as the annual basal area 


increment, at breastrheight say. 


We integrate this eguation by separating the variables— 


da 


1.€. BA — fa — di 
H oe Er eie 
Ë 1 d(BA ja) 
— kl (6 
7e. F; j) SEER H 
16 E j) loge PA —fa)— ku (1) 


Hiati-—ia— dywe have 
D loge BA — fa) — Ry HC 
1 


F loge (BA — fa) — Kt 


1.6. @ — EE 


Substituting in (1) we have— 


A 
1 


— log BA — fa) — kE— ty) ; loge (BA — fag) 
ba N R 
i.e. loge TER EE fe) EE ef G ”) 
—Af A—E) BA — fa). 
ees ENT R 
From this we get— 
rare 
Substituting in (A), or just differentiating (B), we get— 
s ai RGA fa EA EET 


! . 1 : 
Now the growing-space A is egual to ; where n is the 
number of trees per acre. Hence we have— 


ED 


If the growingrspace is expressed in similar units to the basal 
area, which is usually in sguare inches— 
6 -27964 xx 106 


A — sg. inches. 


Finally we thus have for the annual basal area increment in 
sguare inches— 


— (6 D7264 10% — faam) ell SO (CO) 
sd ES! da. da, 
if at; — is we have BA Er 
: da, h 
ss EE -D oe 
Then (6. 6-27 X 10%. fao) 
da da, — k—t) 
Hence Zi — TT Fi ae s ( ee Fo) id (D) 


Ë 


da, . ` : 
Here 4 is the annual increment put on the year just 
previous to suppression and is thus of fundamental importance 


in determining the future basal area increments of the trees. 


Integrating (D) again we fnd— 


— Di N. dag. —BE—L 
say EE (1—e ) Map (EE) 


This is a more convenient form than (B) as it contains terms 
which are more directly measurable and admit of a direct 
physical interpretation as will be shown later. 


Eguation (C) gives the relation which the annual increment 
da. 
Ge in any stand bears to the following — 


(1) The density of stocking. 


(2) The soil guality and tree species as expressed by kb 
and kf. 


(3) The time, provided : 2 &. 


APPLICATIONS OF THE DIFFERENTIAL EOUATION 
OEWEROMT EI 


A. The Free-growing 'Tree. 


For a freegrowing tree A is not limited but increases 
with the age of the tree. Hence DA increases as fa increases 
but we may take it that as soon as the tree attains its full 
vigour the difference (bA—a) remains constant. Referring 
to the differential eguation (A) this would mean that after 
a time— 


da 


ER BORSEL — GM SERP og soe PD PA E (1) 
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"This is in accordance with the investigations of R. Weber (6), 
who found that after an initial period of slow growth, the 
sectional areas of the annual rings of a tree, whose growth 
is unchecked, tend to remain the same up to an advanced age. 


Integrating eguation (1) above we get— 


2 se mt Ee Me Ke N EE N, 
ifat # 


Si 4 — da then 
AG — mi HE 
Cc — dy Mi 


Substituting for c in eguation (2)— 


As mt aa BA RE 
where a — Basal area of the tree at B.H. say 
# MA ge oi iheltreë 
 — Initial period of slow growth 
aA — The total basal area of the initial growth period at B.H. 
m — Constant annual basal area increment. 


lt is dificult to obtain a tree in practice, which is absolutely 
freegrowing, unless a small stand is planted at a wide espacement 
with this specific object in view. As such a stand is not 
available at present, the authors examined an individual tree 
growing more or less alone. 


Table 1 gives the data obtained from the stem analysis 
at B.H. of what the authors consider a free-growing P. pinaster 
tree, growing on poor soil in Stellenbosch, C.P. lt will be 
observed that there is an initial period of slow growth up to 
the sixth year and fairly constant basal area increments from 
the seventh to the thirtieth year. 


The mean basal area increments for the 24 years 
— 4.687 sg. imnches. 
Standard deviation - .45 


Hence 68% of the annual basal area increments fall within 
4-6874-45 sg. inches, ie., a variation of 4 96% from the 
mean. Schumacher and Meyer (5) found in a study of 12 
dominant white fir trees that the #luctuation in basal area 
increment, due to cdlimatic influences over a period of 63 could 
be expected to be 9:5%. 


MA ANBIEE ST 
Age in Basal area sg. ins. | Annual increment, 
years. Measured with sg. ins. 
Planimeter. 
1 155 

2 “503 -348 
3 1-852 -7A9 
i 1:266 914 
5 4360 2-094 
6 7:-995 3-665 
7 12-211 4.286 
8 17-198 4.987 
9 22-600 5 -402 
10 27-371 4.771 
1 382-214 4.843 
12 436-943 4.799 
18 41457 4.514 
14 45 -871 4.414 
15 50757 4.875 
16 54-814 4-057 
17 58-814 4-000 
18 63-014 4 -200 
19 67 -943 4.999 
20 71-957 4-015 
21 76 -700 A.743 
22 80 757 4 -057 
23 85 -028 4.971 
24 89 - 885 4.857 
25 94 .976 5-086 
26 100-071 5-100 
27 105-400 5-399 
28 110-100 4.700 
29 114-300 4-200 
30 119-414 5-114 


APPLICATION TO UNIFORM EVENAGED STANDS. 


The above formula applies to a single tree, the growth 
of which is not limited by insufficient growing space. If we 
want the law of constant basal area increment to apply to a 
uniform evenaaged plantation, the trees constituting the stand 
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must be freegrowing. ln practice this is defined by O'Connor 
as follows — 


“For practical purposes, trees in uniform evenraged stands 
of different densities in the same locality are regarded as free 
growing, while there is no differentiation in the mean trees of 
the different stands. 'The moment the mean tree in a stand 
of a particular density starts to exhibit diferentiation in crown 
or stem as compared with the mean tree in a stand of lower 
density, then the trees in the stand in guestion cease to be 
freegrowing and are starting to compete with one another.” 


It is therefore possible to obtain in practice, stands of 
free-growing trees by means of early freguent and judicious 
thinnings. 'The above formula will then apply to the mean 
tree in such a stand. 


The growth of evenaged coniferous plantations with the 
growing space of the trees restricted will be examined in the 
second part of this paper. 


B. 'Trees with Restricted Growing Space. 


In deriving the differential eguation we took it that due 
to mutual competition there was a deficiency in annual food- 
supply which got more intensive with time. 'The annual 
food-supply being constant, viz., bA, but the demands of the 
growing tree increasing, the demand obviously must exceed the 
supply at some time. 'This time we can identify with & in 
the eguation (C). lt is the time at which the limitation in 
growing-space makes itself felt. 


THE RELATION BETWEEN THE PASAT SAREA 

INCREMENTS IN TWO SUPPRESSED EVENAGED 

STANDS OF DIFFERENT DENSITIES AT ANY AGE 
COMMON TO BOTH. 


In two stands of densities m and n: let (), and (&)s be 
the time suppression starts; then at any future age t common 
to both stands— 


k 


(2) NE OE EE OE OE 


' da 
EE) ns (6:87 N 10% b — ag, me Von) 
(2) m (6:27 X 168 X b— las): He EE 


(6:27 X 10% X 6 — fam) RUE — os) 
MEESTE T ERK 


It will be shown later, when discussing the initial period 
of growth, that (“hm — (“oja?%s. This fact makes the two 
expressions in brackets the same, and we thus have— 


EE aie GR) — (toe) 
EET. 
Ed re ii where 4 — (Ak — ole 
Hd 


In stands of the usual planting distancest (lok — (oa) varies 
from a few months to a few years. It will be shown later that kf 
is a small fraction of the order 107. If x is also a fraction, ie. 
less than a year, then— 


— va Me 
EE AE N 


Hence for two even-aged stands of densities n, and n: in which 
Suppression starts during the same year— 


(se 
Abe ers tor practical purposes. 


sa) My 

EG). 

If suppression started one or two years later in stand nm, than 
in n: Where n: 7 n— 
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G) 
GEE ye 
the ds ei os and “e- ae respectively. 


(ES) ed n 
diJ, 


It can therefore be considered as a fundamental law that 
in two even-aged stands of different densities, both under 
suppression, the basal area increments at any subseguent age, 
common to both stands, will be directly proportional to the 
mean growing spaces of the trees in each stand, multiplied by 
a correction factor. 'This factor arises owing to the difference 
in the times taken by the limitation in growing space to make 
itself felt in each stand. Where the difference is less than one 
year the correction factor can for practical purposes be omitted 
or rather taken as unity. 'The growing space is expressed as 
the reciprocal of the number of trees per acre. 


This law may be stated differently: if we take the age of 
the two stands to begin from the moment suppression starts 
then the annual increments in the two stands at the same new 
suppressed age are in the ratio of their respective growing spaces 
Ee the correction factor. 'This will be shown graphically 
ater on. 


In practice the simplest application of the law will be to 
two stands of the same species, planted at different densities 
and then left unthinned. 'The law will apply to the increments 
in these stands after suppression has started. Since plots have 
never been laid out with this particular object in view there 
is no data available for the straight-forward testing of the 
law. 'There is, however, data available from plots of different 
densities established by means of thinnings and the authors have 
perforce to avail themselves of these. 


Now, if a stand is thinned, then the trees are freegrowing 
for a period, relative to their condition just prior to thinning, 
as will be shown later on. If left unthinned the trees soon 
become suppressed again. ln this new suppressed state the 
increments may be compared with those of any other suppressed 
similar stand of different stocking. 'These circumstances obtain 
in the frst example which will be used to demonstrate the 
validity of the law. 


ld 


Example 1. 


Use will be made of the data in Table II which are 
abstracted from data of sample plots of P. patula at Weza 
Plantation, Natal, given by O'Conner (4). 


TADEEAIIS 
Plot 5 Plot 2 Plot 4 
Density 959 343 183 


M. Diameter at B.H. 
ins. AS FM] 


M Baal area BE. 
EE. N oSzS Midsd N is6ed || Pied || A1364 Ma7 


dnetemengt As Mer -0981 -0800 1012 


The plots received the following treatments— 


Plot Original Stock. Treatment. Stock at 
No. 63, years. 
5 952 trees per acre No thinning 952 
2 leid dd Thinned at 86 years to 343 


343 trees per acre 


4 OM EE, Thinned at 38 years to 183 
343 trees per acre, and 
again at 5 years to 
183 trees per acre 


ii 


The reduction of the original stock in plots 2 and 4 having 
been done to avoid competition, it follows that the stands were 
already competing at 3 8/12 years. Let us assume that com 
petition started at 3 years. 'Trees in plots No. 2 and No. 4 
went through what may be termed a period of response, or 
relative free-growth, but according to O'Conner (4) competition 
started again at 5 5/12 years when plot No. 4 was further 
thinned. For the period from 5 5/12 years to 63/12 years 
plot No. 2 is thus under suppression and the basal area 
increment should therefore fall within the scope of the eguation. 
In other words it is assumed that after the period of response 
the trees in the thinned stand are in the same condition as 
trees would be in a stand, the same age but started originally 
at the espacement of the thinned stand. If this is the case, 
as will also be shown theoretically later on, then the increments 
from the 5th to 6th year in plots 5 and 2 should be directly 
proportional to their respective growing spaces with the 
correction factor. | 


kf is not known for this stand but, as has been stated 
previously, is of the order 10” for any species. Xx in these two 
stands must be less than 5 5/12—3 8/12, ie., 18 years, possibly 
less than 1 year, and in that case can be taken of the order 
unity. 


Hy -O9S1 343 


— 


ë ( 
Thu ME dv se. 
N TED ia 08 — 958 


or. :a5il *aoo ME EO 


This gives a variation of 2.5% from the result we expect 
and can be considered as near enough to be taken as proof 
of the validity of the application of the law to these two stands. 


N.B.—The trees in plot No. 4 are freesgrowing and hence 
Outside the scope of the formulae except where the stocking 
consists of the maximum number of freegrowing trees the 
area can sustain at the particular age under consideration. 'The 
law can therefore be applied to calculate approximately the 
maximum number of free-growing trees possible on plot No. 4 
at 6 2/ 12svears 
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d 
e.g. ( si increment of plot No. 2 with density ?%, 
1 


No. 4 ,, s Ns 


GO... 
AE 


343 X -0S 
“1012 
271 trees approximately. 


I 


Taking plot No. 5 and plot No. 4— 


952 X -0281 
1012 
— 264 trees approximately. 


9 EE 


The maximum number of freegrowing P. patula trees at 6 
years on this particular site will therefore be 4 270. 


Example 2. 
dADEES HI 
MS PANGS SP No PiNos 
374 trees/acre | 496 trees/acre | 258 trees/acre 
1928 | 1936 |. 1928, | 1986 | 1928 | 1936 
Mean B.A. sg. it. ... | “094 “094 “043 “076. | “0517 || -099 
Inererment Aa “040 “03% -048 


The data in table JIT, see Chandler (1), gives the mean 
basal area at BH. of Eucalyptus Astringens in sample plots, 
laid out on an area of natural regeneration containing 2,089 
trees per acre on reserves near Wickepin, Western Australia. 
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From the nature of the establishment of the forest, it must 
comprise some trees in which limitation in growing space made 
itself felt sooner and others in which it made itself felt later 
due to the natural unevenness in stocking. 


In thinning the trees on the plots in 1928, however, the 
trees were selected in such a manner, that for those retained, 
due consideration was given to the general vigour and espacement 
and those removed included the suppressed and dominated. 'The 
suppressed and dominated trees were most likely those in which 
suppression started first and the vigorous growing ones those in 
which it started later. Hence due to the nature of the thinnings, 
those left would be mostly trees of the same class and we 
could assume that for the trees comprising the plots, limitation 
in growingspace made itself felt more or less at the same time 
and we can again take the correction factor as practically unity. 
Furthermore, if we compare plots No. 1 and No. 2 for the 
incremental period 1928-26 exactly similar relative conditions 
as regards period of response and suppression obtained and 
their increments for this period should therefore be directly 
proportional to their respective growing spaces. For nm, and ns 
the number of trees at the end of the period was taken. 


7), 
n) 


je dd EE IEGOS 


— dé becomes Ed — E 
1 


-(032 374 


Now ( 
IK 


ie., a deviation of 9.6%. 


In plot No. 3 the trees were free-growing for the incremental 
period. Using the values of the mean increments and respective 
densities in plots No. 1 and No. 2 and the mean increment in 
plot No. 3 the approximate number of freegrowing Euc. 
astringens trees the area can support, at the particular age of 
the stand, can be calculated in the same way as in the previous 
example. 
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THE RELATION BETWEEN SUCCESSIVE BASAL 
AREA INCREMENTS IN A SUPPRESSED 
EVENAGED STAND. 

Let t. be the time suppression starts in an evenraged stand 
with density of stocking n which has been left unthinned: 
Ee P subseguent successive basal area increments will be given 

VY 


da k —AA— 
es (b X 2-7 X 10% — fnas)e TE) 
da k —RIE—L 
Er se (bx 9-7 Xx 105 — fnas) e Te —ty) 
da 
EE MAR) 
day —RE—L) 
dié 
kf 
sé where ($—t) — one year 
da, 
sm di Ef 
S s 
oimilarly da. ë 
da, 
di kf 
WE aM ER 
di 
da, da; das da, 
Es DEP EE TE dt 
dk n OE REBEL LAS: RE it 
di di di dt 


If our assumptions are correct, then in an evenvaged stand 
under suppression, if left unthinned, the basal area increments 
in successive years should decrease according to a geometric 
progression. 

To obtain factual proof of the above law three P. pinaster 
trees were cut in a 27 year old stand in Union Park, Stellen: 
bosch, C.P. 'They were planted at an espacement of 9 ft. by 
9 ft. and left unthinned. 'The trees were selected so that two 
were about the mean D.B.H. of the stand and were codominant, 
and one a large dominant tree. 

The data obtained from the stem analysis of these trees 
are given in table IV. From the measured mean diameter in 
inches of the successive annual rings at B.H. a smooth curve 
was plotted against the time in years for each tree. 
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See graph No. I. From these curves the mean diameters 
were read and used for calculating the basal areas. 'The basal 
areas thus obtained for each tree were again plotted against 
the time. See graph No. II. From these curves again the 
basal area increments were read. 


DIAMETER IN INCHES. 
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Basal area CcurTves. 


Our formula for the basal area increments is— 


—y) 


od N - 
If we take logarithms on both sides— 


— loge (GO 9:7 105 — fa) — RIT 


da 
ER 
This is the eguation of a straight line if we plot log: se against t. 
Graphs were constructed accordingly by taking the Napierian 


IM; 


logarithms of the annual increments as read from graph TT, and 
plotting these against the time. See graph II. 


Ge 
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Graphs of the logarithms of the increments against the time. 


These graphs supply clear proof that from the moment 
suppression sets in the logarithms of the increments plotted 
against the time give a straight line relationship. 


It can therefore be stated as a fundamental law that in 
an even-aged stand under suppression, if left unthinned, the 
basal area increments decrease according to a geometric 
progression. 
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We may therefore regard our formula as proved in so 
far as the results predicted conform to actual observation. 


INTERPRETATION OF GRAPHS. 


A. 'The Initial Period. 


With reference to graph No. TT it will be noticed that 
the initial freegrowing period which falls outside the scope of 
eguation (D) is also represented by straight lines, with an upward 
slope, however. 


To account for this feature we shall have to derive a 
growth-formula for this initial period also. 'This period differs 
from the period of constant growth of a mature freegrowing 
tree, for which we have already accounted. Jf we assume that 
for this initial period the basal area increment in unit time 
is proportional to the basal area at any given moment we 
will have that the growth da in a time dt is given by— 


da ei eadk ae ET EE 
where c is a constant of proportionality. Integrating we get— 


loga—sd-d-d 


where d is a constant of integration. 


diatis is 


loge a, — dd 
a 
or Me VY dor SE N 
da c(i— 
Hence mae ë ho). Ms Re EA 6 
aa da, clt—t) ) 
or sé Mm) 
. da 
. loge TE logease-kelisi) Es, 


d | 
If we now plot log EE against t we get a straight line, as is 
actually obtained in our graph. 


id 


Each tree in a stand will have a graph of its own which 
need not coincide with that of another tree in the same stand: 
it will, however, have the same slope since we take the constant 
Cc to be the same for all trees of the same species in the same 
locality. 


This means for a given species in all stands in the same 
locality irrespective of their degree of stocking, the logarithms 
of the annual increments of the mean trees, during the initial 
free-growing stage, when plotted against the time, will give 
one straight line of a given slope and position. 'The length 
of the line for a particular stand will be determined by the 
duration of the initial period of free growth, which again is 
determined by the density of stocking. 'The smaller the density 
of stocking the longer the line will be until eventually it bends 
over and becomes horizontal for a stand which is free-growing 
for the entire life of the trees. When the density of stocking 
is such that suppression starts during the initial period, the line 
is only defined in terms of our edguation (5) up to the time 
suppression starts. When suppression starts it means that the 
entire ground is occupied, i.e., the productive capacity is tested 
to its utmost limit. In the last year just before suppression 
starts the total increment in a stand of density n, is thus the 
same as in a stand of density ns, ie.— 


AE ETED] wo aa TE AG 
but ( )- c(a), and ( )- c(a). trom (]) 


These eguations hold for the whole of the initial period, 
hence also for the last increments— 


je. ( ie ) c(ay) and ( “re le oa) 


n ( oi ' — Mm Elaoh and % ( 7 ). TIGE 


di 


2) 
Substituting in (6) we have 
n (ak — Ms Clay) 


ie. Mm (ah — Mg (ay) 
This proves our former assumption on page 9. 
B. The Period of Suppression. 


Referring to our eguation (D) we have on taking logarithms 
that 


da da 
loge — — loge Ee ESME EE BEE Ed 


This is the eguation of a straight line if logs da is plotted 


against t and its slope is kf which is an absolute constant for 
a particular species in a given locality. 'This means that for 
a given species in all stands, under suppression, in the same 
locality the logarithms of the annual increments of the mean 
trees plotted against the time will give a series of parallel straight 
lines one for each degree of stocking. 


FIGURE 1 


In figure 1 above is illustrated the history of the growth 
of the mean trees in four even-aged stands of different degrees 
of stocking ni, Ns, ns and n. per acre. 'The initial periods of 
free growth all lie along the same line AB and the subseguent 
successive periods of suppression branch off from this line at 
the points a, b, c and d, which signify the times at which 
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suppression started in these stands respectively. In the limiting 
case, with such a density of stocking that the trees are free 
growing for their entire life we will have the line ABC where 
the latter part BC is horizontal and represents the period of 
constant basal area increment. 


The numerical value of kf is obtained as follows: In any 
given stand— 


loge oe — loge Ee EE —) 

log se ore BE —) 

loge s loge — BE—R) 
 loge EE loge 

T er. 


Referring to graph No. IIT and substituting the values 
obtained from the graph in the above eguation, we get the 
following table — 


da, da, hi 
free é is log FT log FT ke] | 
A 7 12 1 “105 “360 -143 
B 7 12 1 -980 “205 155 
é 7 12 1:353 660 -135 


Average &f — “144 


We may now also calculate the value of kb which will be a 
site corefficient. 


We have— 


da od oe 
di 7) 


B X 6:2 X 106 — fame HY) 


"Taking tree A which lies between trees B and G to represent 
Mk d 
the mean tree of the stand more or less we take its ed at say 
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the 14th year and its as at te, ie, 6.2 years. From the graph 
we can find the values and on substitution we have— 


Eb —kf( 
` —— EE -9 D k (eg) 
1-116 sig SAT 10 — Af )e 


14—6:9) 


7:8 kf 
Bb (1-116e - Af X 12-3) 540 
Ë 6 (OS 
Now Af — -148 for tree A 


Finally on substitution of this value for kf— 


kb — -000465 


Thus s — 310 and o — -00323 


C. The Period of Release. 


Referring to graph No. II we notice a sharp. change in 
the slopes of all the lines at some time during the suppression 
period. 'This change is especially marked in the case of the 
dominant tree C. 'The upward sweep in the line of tree C 
is so marked as to make the new line apparently parallel to 
the initial period of free growth. We may therefore suspect 
that this period corresponds to a new period of free growth 
due to a sudden change in n. If this is the case an eguation 
similar to that of the initial period will apply. 'Tree C as 
previously stated was a dominant tree. lt was surrounded in 
the stand by several badly suppressed trees and some trees 
bordering on it had actually been removed a few years ago. 
This release was so marked as to arrest suppression and render 
the tree relatively freegrowing again, hence the upward turn 
of the line. 


'To derive the differential eguation of this growth period cd 
on the graph we proceed as follows— 


The eguation of the line cd is obviously 


da da 
loge Pe loge 4 elf—me ak N N iN 
é : da da, clt—h) 
From this we get or Ee EE ë of KM ve NE 


Integrating we get: 
EI 6 Hd oe EN ii) 


TÉ at PES Fe d — My 

dl da ed) 
EN — Erg 

C 


th — oe 
En G d 
N el ER 
Er AM EE. DR ES 
1 da ct—h) 1 da 
H ei RS 1 1 
ence a EE Ha, EE (5) 
Substituting from (2) we get; 
1 da da 
ME id Ed. ad. 1 
Ee MR ie Hi 
by da # da, 
di EE 
Substituting from (4) for a, 
da 
— —— dl) ) 
FT c(a Dier vars “ri Ge ant (6) 


This differential eguation describes the growth for the period 
of release from t, onwards and takes the place of 


iE — ca for the initial period along AB. 


SE | 
“Toe $nd a and hence GE in terms of the time we have to 
integrate (6)— 


da 
EE ELE E— 4 
J a—d f AE 
's lege ad. d HE 


Ti atist dd 
loge (a —ddAE 


On substitution for E 
B sd 

e — —l 

log GE Es 5) AG N) 


ora—dsl@A—d)e 


From this 1 — ca — dye (8) 
AE (5) 
EE 
where Ee NE the value of at VR 
dié dié 


'Taking logarithms on both sides of (9) we get our original 
eguation (1) which gives us the line cd. 


If we compare (1) with the eguation for ab, VE— 


da da 
loge eri. loge GE H CE—ty) 
da, da, f 
and take loge EE loge Er and also c(—ty) — ct—&) 


(where the ts are naturally not the same) we see that the 


log E `s may become egual. In other words we have the 


remarkable fact, which should also be evident by a study of 
the graph, that the free-growing tree in the stage cd, although 
much older and thicker than the same tree in the free-growing 
stage ab, actually puts on the same increment per year as it 
did at a certain stage of its life in the period ab. 'The years 
of suppression have actually had the effect of putting it back, 
as far as growth is concerned, to a period in its initial free 
growing stage arrived at by drawing a line parallel to the 
traxXIs. 


The line cd appears to be parallel to the line AB and if 
this is the case where complete release is affected through 
thinnings it would mean that the logarithmic increment in the 
rejuvenated or released stand is the same as for the young tree. 
In case the line cd should be found not parallel to ab it will 
not change the form of (6) nor will it affect our subseguent 
development in its essentials. 


We must now settle the guestion: How far up will the 
line cd go? Tt will turn sharply down the moment suppression 
sets in again. When will suppression set in? 
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When the total annual basal area increment has reached 
the same amount as would have been attained in the final year 
before suppression by an unthinned stand of the same density. 


We have no data at our disposal to prove this contention, 
but suppression in the released stand can only set in when the 
trees occupy the entire area again and limitation in growing 
space makes itself felt; at this moment the increment in the 
thinned stand must be the same as that in the year suppression 
starts in a stand which was originally planted at the espacement 
of the thinned stand. 


With the line representing tree A the upward sweep is 
slight compared with that of C. 'This would correspond to 
a partial release of suppression obtained by suppressing its 
neighbours. 'The release was so slight that the suppression set 
in during the same growing season again and the line actually 
proceeds in small zigzags. lt should mean that in course of 
time tree A, which is now still codominant, would become a 
dominant tree if the stand is left unthinned. 

With the line representing tree B a .downward sweep is 
encountered. 'This would correspond to an increased degree 
of suppression and would mean that in course of time tree B, 
which is now still co-dominant, will become entirely suppressed 
if the stand is left unthinned. 


The graphs of A and B indicate that if the stand is left 
to nature, tree A will survive and make up one of the trees 
in the stand at maturity, while tree B will become more and 
more suppressed with time and might eventually die of. 'The 
tree is loosing its vigour resulting in an increasing rate of 
decrement and from this stage falls outside the scope of our 
formulae. 


THINNING SCHEDULES. 


It will now be shown how practical use can be made of 
the graphs of log: ) against time as a guide for drawing up 
thinning schedules. 
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da, 
Sr ( di 


with density %,. 


then sa) ad 1] (6 2 Da we E — fa ee N 
1 dl 


di 
Let ( 


with density 7. 


) be the increment one year after Suppression in a stand 
1 


da i Ek ks 
- 7) be the increment one year after SUppression in a stand 
& 2 


1 al 
da k —HWE— ) 
er (Ge EE. (BOX 2-7 105 — Fa) mg)e  eG 


(2) Th (gl 1 
di J, as EE RUE —ty) EE) 
( da, ) Mm 
diJ, 
Put ed E—) —1 


GP on 


and 


€ EK 5 a) ds é 
or log (G ) log ( ia) oe loge n loge , 
' GRAPH IE. 


KI 


Mg 
FREE GROWING 
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Thinning schedule graphically illustrated. 
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If sufficient data is therefore available to draw a graph of 
log: against time for any one stand of known density, then 
graphs for all other stands of the same species in the same 
locality and on similar soil can be constructed. 


The graph representing tree A in graph No. III is taken 
to represent the growth of the mean tree of the stand investigated. 
The density of the stand is known to be 540 trees per acre and 
hence all the data is available to construct a series of lines 
representing stands of P. pinaster of different densities growing 
under similar conditions. 


These are represented in graph IV which is constructed 
as follows — 


The graph representing tree A in graph No. TT is trans 
ferred to graph No. IV and represented by ABC. loge one 
year after suppression, i.e, in the 7th year, is read from the 
graph at point G. lt is egual to 1-12. 


The line DGE is drawn through G parallel to AB. 


Let us now calculate the l1st increment after suppression 
for a stand with 400 trees per acre— 


loge (ze). — loge (sa ).- loge n; — loge 


' da, N 19 
*. loge ( ). — loge 540 — loge 400 4 1-12 
mr TT 


The point L representing 1-4201 is read off on HG produced. 
A horizontal line is now drawn through L to cut DGE at K. 
The line KN drawn parallel to BC represents the line for a stand 
with density of stocking 400 trees per acre which has never 
been thinned. 'The table below gives calculated logarithms of 
the first increments one year after suppression of a series of 
stands of the densities indicated. 


Density | 1000 s00 | 600 400 300 200 106 


da 


FT -3038 | -7970 | 1-0146 | 1-4201 | 1-7078 | 2-1143 | 2-8064 


loge 


The lines representing stands of the densities indicated above 
are constructed in the same way as that for density 400. 


Now the average maximum basal area increment a free 
growing P. pinaster tree can put on in the locality was shown 
to be 4.67 sd. inches, see Table 1. 


In practice therefore the graph representing a stand with 
the trees freegrowing for the entire rotation will follow the 
curve ARS. 'The graphs for all stands with densities higher 


than that of the freegrowing stand must fall below the curve 
ARS. 


Theoretically some of the graphs go above the line RS, but 
in practice due to the limitation of growth of a tree at the 
maximum for the particular site and climatic conditions the 
actual points obtained by observation must all be below the 
line RS. 'The theoretical portion above RS, however, enables 
us to get the true position of the actual lines representing the 
stands of the different densities indicated. 


Take as example a stand with density of stocking 200 
per acre. 'This stand will theoretically follow the line APM, 
but if actual data could be obtained from such a stand the points 
will follow the line ARTM, i.e, the stand will be freegrowing 
for the entire initial period but its growth along AF will be 
arrested, as soon as the maximum growth is attained at R. 'The 
stand will now have a period RT of maximum free growth 


after which suppression sets in, and if left unthinned, it will 
follow the line TM. 


The series of curves constructed may be looked upon as 
derived from fictitious stands of the different densities. 'These 
shadow stands, so to speak, enable us to arrive at values for 
real stands of given densities which have been obtained from 
a single stand by successive thinnings. 
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Let us now consider an evenaged stand planted at a 
density of 1,000 trees per acre and work out a 5 yearly thinning 
schedule. 'The initial period of growth of this stand will follow 
the line AF and suppression will start at the point “a” corres 
ponding to an age at BH of 3.6 years. During the early years 
of the stand a certain amount of suppression is considered 
beneficial, both from the silvicultural point of view and the 
point of view of the guality of the timber produced. 'The stand 
is now allowed to be suppressed up to the 6th year, i.e. 10 years 
from date of planting. At this stage the stand has reached 
point b and is then thinned to 800 trees. lTmmediate response 
will now take place along the line bc, as has been shown before. 
Response stops at c where c is the same height as point x of 
the shadow stand for 800 trees. 


Suppression now sets in along line cd and if we now thin 
to 600 trees per acre at 11 years BH age, ie. 5 years after 
the previous thinning, release sets in again at d. 


This zigrzag process is repeated, the thinnings taking place 
at 5 year intervals, 200 trees being removed at each thinning. 
At the 21st year we remove only 100 trees; the release takes 
place along the line hi where i is on the maximum freegrowing 
line. - We, however, project the line to the point j, which is 
on the same horizontal line as the point w of the shadow stand 
representing 300 trees. From j we draw a line parallel to all 
the previous suppression lines to where it cuts the horizontal 
free growth line at k and from there to the point | which is 
on the 26th ordinate. 'The curve of growth we assume to 
be the line hikl. 'This is an extra-polation which in actual 
practice may perhaps not be borne out as the suppression may 
start sooner at some time less than k but greater than j. In 
this case 1 will just be displaced somewhat lower, whereby In 
is displaced to the right, while o may hardly be shifted at all. 
In the end our main result may scarcely be affected by this greater 
finesse. (Only actual experiment over a period of years can 
decide this point. 


At 1a further 100 trees are removed and the curve now 
followed is given by lmo; at o, which in this case is during the 
34th year, the stand is thinned to a final density of 100 trees 
per acre, which will start suppression again at g at the 49th 
year. 


ap 
It will be noticed that the stand was freegrowing except 
for the short period klm, from i to g., ie, from the 22nd to 


the 49th year. If the final Crop is now removed we will have 
a rotation of approximately 55 years from the date of planting. 


TABLE V. 
THINNING SCHEDULE 


Initial stand 1000 trees per acre. 


No. of Basal No. of Basal Diameter 
Age from trees re- area, trees area BR. 
Planting | moved in removed left in left sd. of mean 
thinnings | sa. inches stand inches tree 
10 200 brig) BOO 6920 2.40 
Es 200 3636 6OO 10908 4.90 
20 200 6190 400 12380 6.29 
25 100 500 300 15000 7.98 
20 100 7300 200 14600 9.64 
238 100 EDE 100 11272 11.98 
35 — —— 100 18044 15.14 


Table V, which is based on graph IV, forecasts a complete 
record of the trees removed during thinnings with the diameter 
at B.H. of the mean trees. 'The total basal area removed in 
thinnings amounts to 35,168 sg. inches. 


The final stand will have a total basal area of 18,044 
sg. inches per acre, with the B.H. diameter of the mean tree 
eguivalent to 15.14 inches. 'The above figures are for a stand 
On very poor sandy soil with a rainfall of 25—30 inches. 'To 
grow P. pinaster trees to a larger size on this site would mean 
a considerable lengthening of the rotation by a further thinning 
at 55 years if a new period of suppression is to be prevented. 


The thick line in graph IV represents the growth line of 
the mean tree in the thinned stand. It is put forward tentatively, 
largely on theoretical grounds. 'The thinning 'schedule is 
hypothetical and just an illustration of how our eguations can 
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be applied in practice to draw up thinning schedules for certain 
areas and tree species, with a minimum of sample plots, from 
which the reguired data can be derived. Should our theory 
be verified in all its conseguences by actual experiment, it will 
place the art of thinning in evenraged stands on a definite 
scientific basis, where guantitative prediction will take the place 
of empirical estimate. 


SUMMARY. 


1. A differential eguation of growth is derived giving 
the relation between increment, growing space and time for 
a given species in a given locality. 


2. As a limiting case this eguation is applied to the free 
growing tree. 


3. It is established that for two suppressed evenaged 
stands of different densities the basal area increments taken at 
the same age from the date of suppression in each stand are 
directly proportional to their respective growing spaces. 


4. 'The increments in an unthinned even-aged stand under 

suppression decrease with time in a geometrical progression, i.e., 
a se 

there is a straight line relationship between log. 3, and the time. 


5. With the aid of data obtained from a single local stand 
graphs are drawn for a series of stands of different densities 
giving the relation between logarithmic increment and the time. 


6. With the aid of these graphs again a new graph of 
logarithmic increment is obtained for a stand which is succes” 
sively thinned out at definite intervals. 


7. From this last graph a tentative thinning schedule is 
drawn up. 


Finally the authors wish to thank the following bodies and 
persons for facilities placed at their disposal or helptul criticism: 
The Stellenbosch Municipality, Mrs. J. H. Marais of Coetzen 
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burg, Mr. Cornish-Bowden of Lourensford, Somerset West, and 
Prof. E. ]. Neethling of the Stellenbosch-Elsenburg Agricultural 
College. 


ADDENDUM. 


The above paper was handed in for publication in December, 
1937. In May, 1938, a dissertation by Walter Peschel reached 
the one of us (A. C. Cilliers), in which a critical examination 
of all growth formula's for trees was given up to date. 'The 
dissertation was published in the `“Tharandter Forstliches 
Jahrbuch,” Hetft, 3/4, 1938. 


Peschel derives a growth formula which has been tested 
for the height growth of trees, which superficially agrees very 
closely with our own for the basal area growth of trees. 


He assumed that the increment in growth is proportional 
to the difference between the maximum and present value of 
the growth-variable— 

dy 


re bk REEM N IR EN 


From this he gets by integration that— 


de EE C i 7) EE Ns 
€ 
This may be written— 
"ER VA SE ER EN 
We had for the basal area— 
aratg TE Ure MEN ie EE 
san de AE RE 


The similarity between (C) and (E) is obvious. He uses Xx 
for the time, where we use t. 
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Our own formula was, however, derived in a totally 
different way, and, to my mind, in a physically “einwandfreier” 
way. lt is, moreover, much more general in that it takes into 
account the fact of suppression in its dependence on density of 
stocking. 


Peschel's annual increment would be— 


and our's would be— 


da k 


is BEL) 
EK (67264 x 10% — fa, me Me 


where n is the density of stocking. (Our constants kb, kf, as, 
te are easily determined after about five years of growth, 
whereas the Ymax of Peschel is a value only determinable after 
a time comparable with the life of the tree. With us kb and 
kf are soil and species constants, easily derived by a stem 
analysis. 


Finally, our formula enables us to draw up a (tentative, 
pending further experimental proof) thinningsschedule, which is 
the real purpose behind the search for a growth formula. 


In so far as it goes, however, Peschel's work confrms 
our own. ln his final paragraph Peschel uses these words: 
“Wenn man sich eingehend mit Wachstumsvorgangen ver 
schiedenster Art beschaftigt, so muss man zu der Ueberzeugung 
kommen, dass ein allgemeingiltiges Gesetz waltet und dass es 
nur eine Frage der Zeit ist, bis es gelingt, dieses Gesetz in 
eine mathematische Form zu bannen.” 


I hope that we have succeeded in Anding this mathematical 
form, or at least a substantial part of it. 


IG EGIELIERS. 


(1) 


(2) 
(3) 


(4) 


(5) 


(6) 
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